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We focus on uncertainties in supernova measurements, in particular of individual magnitudes and 
redshifts, to review to what extent supernovae measurements of the expansion history of the universe 
are likely to allow us to constrain a possibly redshift-dependent equation of state of dark energy, 
w(z). focus in particular on the central question of how well one might rule out the possibility of 
a cosmological constant w = — 1. We argue that it is unlikely that we will be able to significantly 
reduce the uncertainty in the determination of w beyond its present bounds, without significant 
improvements in our ability to measure the cosmic distance scale as a function of redshift. Thus, 
unless the dark energy significantly deviates from w(z) = — 1 at some redshift, very stringent control 
of the statistical and systematic errors will be necessary to have a realistic hope of empirically 
distinguishing exotic dark energy from a cosmological constant. 



I. INTRODUCTION 

Eight years ago two teams observing distant Type la 
supernovae (SNe la) announced evidence that the ex- 
pansion of the universe is speeding up [H, . The dis- 
tant supernovae appear dimmer than they would be in a 
matter-only universe. If this is a true distance effect, it 
implies that about 70% of the energy density of the uni- 
verse reside in a smooth component with negative pres- 
sure, leaving only about 25% is in dark matter and 5% in 
baryonic matter. While compelling evidence for precisely 
this combination was pointed out several years earlier [3| , 
based on measurements of the clustering of galaxies, age 
of the universe, measurements of the baryon and dark 
matter densities, and the Hubble constant, SNe la were 
the "shot heard around the world" as they provided ex- 
plicit evidence for the largest contributor to the energy 
density of the universe. This constituent became known 
as dark energy. 

Since that time the constraints have improved signif- 
icantly: combination of cosmic microwave background 
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experiments [U, large-scale structure surveys [5 
particularly, new SNe la observations ff\. M. M. HQ 
constrain the dark energy equation of state, w 
be within about — 1 ± 0.1. 

Unfortunately, our understanding of dark energy is as 
murky as ever. The simplest model for dark energy is 
provided by the cosmological constant term in Einstein's 
equations, and is still an excellent fit to observations as 
it predicts w — 1 identically and at all times. However, 
even if observations were to pin down the value of w to 
be —1, this gives us essentially no insight into the pos- 
sible source of dark energy. While a cosmological con- 
stant may be the best bet, it is simple to imagine other 
sources of dark energy, including the energy density as- 
sociated with a false vacuum metastable scalar field, that 
would produce a similar value. The only way to get any 
new theoretical handle on dark energy is to be able to 



unambigiously determine a deviation from —1, if such a 
deviation indeed exists. 

Various proposals for dynamical dark energy have been 
put forth which might produce such deviations. Most 
notably, a scalar field rolling down its effective potential 
can provide the necessary ener gy d ensity and accelera- 
tion of the universe 0, III, El, 111 El 111. Generically 
many of these possibilities are already ruled out by the 
data. Among those that are not, none is particularly 
compelling as they also typically do not naturally ad- 
dress the problem why the observed dark energy density 
(~ 10~ 120 Mp L ) is so small, nor why dark energy starts 
to dominate the expansion of the universe only at recent 
times — redshift z < 1. Nevertheless, before we can ad- 
dress such puzzles we need to know empirically if the dark 
energy is measurably distinguishable from a cosmological 
constant. 

Here we will be concerned with two specinc observa- 
tional factors, as well as one overriding theoretical con- 
straint. Observationally we need to determine both the 
magnitude and redshift of individual objects in order to 
map the universe's expansion history. Theoretically we 
have to account for the fact that dark energy has, a priori, 
no predetermined time dependence, and thus our analy- 
ses must allow for arbitrary time variations. 

Essentially all of the consequences of dark energy fol- 
low from its effect on the expansion rate: 
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where m and J7de are the dark matter and dark energy 
density relative to critical, respectively, and we have ig- 
nored the relativistic components. Type la supernovae 
effectively measure the luminosity distance 
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d L (z) = (l + z)r(z) = £ -^- y (2) 

where r(z) is the comoving distance and we have assumed 
a flat universe. Since observations of supernovae allow us, 
at least in principle, to map the expansion history of the 
universe, one can use this data to constrain the nature of 
dark energy. 

This paper is organized as follows. In Sec. [HI we review 
a variety of parametrizations of the expansion history of 
the universe, and therefore ways to measure the prop- 
erties of dark energy. In Sec. IIIII we study the extent 
to which an increased statistical error in supernova dis- 
tances affects determination of the equation of state of 
dark energy and tests of its consistency with the vacuum 
energy value of —1. We conclude in Sec. IIVI 



Many other simple parametrizations of the equation 
of state have been suggested (e.g. [2l|, [22j]); similar pro- 
posals have been extended to the Hubble parameter (e.g. 
|23j). Finally, specific combinations of the expansion his- 
tory parameters, such as the "statefinder" [24j, have been 
proposed as good discriminators between phenomenolog- 
ical dark energy descriptions. It is well worth emphasiz- 
ing that all of the above parametrizations are ad hoc, and 
may lead to biases as the true DE model may not follow 
the form imposed by these functions. Their advantage, 
however, is in simplicity and the fact that two or three 
additional parameters in the dark energy sector will be, 
in the near future, measured to a good accuracy, at least 
when the data from the various cosmological probes is 
combined. 



B. Direct reconstruction 



II. MAPPING THE EXPANSION HISTORY 

A. Key questions and parametrizations 

At the present time, it has become clear that there 
are two major goals that upcoming dark energy probes 
should address: 

1. Is dark energy consistent with the vacuum energy 
scenario - that is, is w(z) = —1? 

2. Is the equation of state w constant in time (or red- 
shift)? 

Violation of either of these two hypotheses would be a 
truly momentous discovery: the former would rule out a 
pure cosmological constant, while the latter would pro- 
vide further evidence for nature of of dark energy via its 
dynamics. 

With this in mind, the most obvious approach is to 
parametrize the equation of state of dark e nergy a s a 
constant piece plus a redshift- varying one [18|,ll9, 20] 
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where the first equation assumes linear evolution with 
redshift, the second is linear with scale factor, and the 
third allows for the transition between two asymptotic 
constant values of the equation of state, with the transi- 
tion at redshift z t with the characteristic width in redshift 
A. Equation^ in particular, has become commonly used 
to plan probes of dark energy as it retains the minimally 
required two parameters and does not diverge at high 
z, while still allowing for the low-redshift dynamics (e.g. 
variation in the value of w(z)). 



Going in the opposite direction, the most general way 
to probe the background evolution of dark energy has 
been proposed in [25|, |2g, |27[ : the equation of the distance 
vs redshift can be inverted to obtain w(z) as a function 
of the first and second derivatives of the (SNe la-inferred, 
for example) comoving distance 
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Similarly, assuming that dark energy is due to a single 
rolling scalar field, one can reconstruct the potential of 
the field exactly 
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where the upper (lower) sign applies if (j> > (< 0) The 
sign is arbitrary, as it can be changed by the field redefi- 
nition, <j) «-> —cj>. 

Direct reconstruction is the most general method for 
inferring the dark energy history, and it is the only ap- 
proach that is truly model-independent (despite some 
claims in literature to the contrary). However, direct re- 
construction comes at a steep price — it calls for taking 
the second derivative of the noisy data. In order to take 
the second derivative, one essentially must fit the lumi- 
nosity distance (or SNa apparent magnitude) data with 
a smooth function — a polynomial, Pade approximant, 
spline with tension etc. Unfortunately, the parametric 
nature of the fitting process introduces systematic biases. 
After valiant attempts to do this using a variety of meth- 
ods for smoothing or fitting the data (e.g. [23|,|28, 29, 30]), 
various authors found that direct reconstruction is sim- 
ply too challenging and not robust even with SNe la data 
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FIG. 1: Example of the direct reconstruction, simulated for future SNa la data and assuming the true equation of state is 
w(z) — — 1. The left, middle and right panels show results when the luminosity distance data are fit with a third, fourth 
and fifth order polynomial in redshift, respectively. Note that, depending on which order polynomial is used to fit the data, 
significant bias, or statistical error, or both are introduced. Adopted from Ref. [28( ] . 



of excellent quality. Figure [T] shows example of the direct 
reconstruction, simulated for future SNa la data and as- 
suming the true equation of state is w{z) = —1 (adopted 
from Ref. [28]). Note that, depending on which order 
polynomial is used to fit the data, significant bias, or 
statistical error, or both are introduced. For an excel- 
lent review of the dark energy reconstruction and related 
issues, see 13111 - 



practice computed directly and here approximated with 
F). Therefore 



F = C 



- 1 = W T AW 
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where the matrix A is diagonal and rows of the decorrela- 
tion matrix W are the eigenvectors ei(z), which define a 
basis in which our parameters are uncorrelated [35] . The 
original function can be expressed as 



C. Principal components 

The next most general method, introduced in [32j (see 
also [33l [3~il]), is to compute the principal components of 
the quantity that we are measuring — the equation of 
state or energy density of dark energy. Principal com- 
ponents are the redshift weights (or window functions) 
of the function in question, and are uncorrelated by con- 
struction. In this scheme, one simply lets data decide 
which weights of the function are measured best, and 
which ones are measured most poorly. The principal 
components form a natural basis that parameterizes the 
measurements of any particular survey. 

To compute the principal components, let us 
parametrize w(z) (same arguments follow for p(z) or 
H(z)) in terms of piecewise constant values Wi (i = 
1, . . . , N), each defined in the redshift range Zj = [(i — 
l)Az, iAz] where Az = z max /N. In the limit of large N 
this recovers the shape of an arbitrary dark energy his- 
tory (in practice, N > 20 is sufficient). We then proceed 
to compute the covariance matrix for the parameters Wi , 
plus any other cosmological parameters such as Qm, then 
marginalize over the latter. We then have the covariance 
matrix, C, for the Wi. 

It is then a simple matter to find a basis in which the 
parameters are uncorrelated; this is achieved by simply 
diagonalizing the inverse covariance matrix (which is in 
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where ej are the "principal components" . Using the or- 
thonormality condition, the coefficients cv^ can be com- 
puted as 

N 

ai = X! W ( Zq ) e i( z a)- (10) 

Diagonal elements of the matrix A, Aj, are the eigen- 
values which determine how well the parameters (in the 

— 1/2 

new basis) can be measured; a {on) = \ i . We have 
ordered the a's so that a(o>\) < cr(«2) < ■ ■ ■ < c(a./v). 

Principal components have one distinct advantage over 
fixed parametrizations. They allow the data to deter- 
mine which weight of the cosmological function w(z) (or 
Poe(z), or H(z)) is best determined. In fact, the PCs 
depend on the cosmological probe, on the specifications 
of the probe (redshift and sky coverage etc) , and (more 
weakly) on the true cosmological model. These depen- 
dencies are features and not bugs, and they make the 
principal components a useful tool in survey design. For 
example, one can design a survey that is most sensi- 
tive to the dark energy equation of state at some spe- 
cific redshift, or study how many independent parameters 
are measure d by any given combination of cosmological 
probes (e.g. [221]). 
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FIG. 2: The four best-determined and two worst-determined 
principal components of w(z) for a future SNe la survey such 
as SNAP [H. Adopted from Ref. JH- 



D. Uncorrelated estimates of DE evolution 

A useful extension of the principal component formal- 
ism is to compute the band powers in redshift of the dark 
energy function, w(z) (or Pde(z) or H(z)). These band 
powers can be made 100% uncorrelated by construction, 
and the price to pay is a small leakage in the sensitiv- 
ity of each band power outside of its redshift range. For 
details on how this is implemented, see [37] . 

The top panel of Figure [3] shows the constraints on the 
four band powers of w(z), adopted from Ref. [37] . assum- 
ing the Riess04 SNa data [3] and a prior on £Im- The 
bottom panel of the Figure shows the same method to 
constrain three band powers, using the newest data SNa 
la data from the Hubble Space Telescope, adopted from 
Ref. [To] , in combination with the baryon acoustic oscil- 
lation measurements @.. The slight preference for w(z) 
increasing with redshift is seen in both panels; however, 
significant systematics still affect the data, in particular 
because of the heterogeneity of the SNa la datasets used. 
Therefore, it is too early to claim any evidence for the de- 
partures from ACDM , and it remains to be seen how the 
results change once we have a systematically more homo- 
geneous and statistically more powerful data set (for the 
requirements on the SNa la systematics, see [38]). 

It interesting to see in Fig. [3] how good the constraints 
are, given that this is current data and that interesting 
constraints are obtained on 3-4 band powers. A superior 
data set with an excellent control of the systematics, such 
as that expected from a dedicated space telescope such 
as SuperNova/ Acceleration Probe (SNAP; [3(|) will sig- 
nificantly improve the constraints, and also allow a finer 
resolution (i.e. more band power parameters) in redshift. 

Nevertheless, it is clear from the existing set that either 
significant redshift evolution in w or significant improve- 
ment in the data will be required before either redshift 
evolution or deviation from —1 could be unambiguously 
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FIG. 3: Top panel: constraints on the four band powers of 
w(z), adopted from Ref. [33], assuming the Riess04 SNa data 
0] and a prior on Qm- Bottom panel: the same method 
applied to constrain three band powers, using the newest data 
SNa la data from the Hubble Space Telescope (adopted from 
Ref. [Ig ]) in combination with the baryon acoustic oscillation 
measurements [f|. 

inferred using this technique. We now turn to considering 
this question in more detail 

III. LAMBDA OR NOT? 

As stressed earlier, one of the most important out- 
standing questions in cosmology is whether dark energy is 
consistent with a cosmological constant. This hypothesis 
can be tested in a variety of ways. The simplest approach 
is to compute a simple likelihood comparison between the 
data and the vacuum energy model w(z) = — 1. A much 
more sophisticated (but admittedly less robust) approach 
would be to perform some type of reconstruction of the 
energy density pde(^) and check whether or not it is con- 
sistent with a constant value (or similarly, to reconstruct 
the equation of state w{z) and check whether it is con- 
sistent with —1). 

Here we consider the A hypothesis test in terms of 
principal components (PCs). We consider models in the 
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wo-w a plane (see Eq. ) and as we outline below, we 
use the best measured PCs to determine if we can statis- 
tically distinguish distinguish the model from ACDM . In 
this section, we assume future SNa la data with 3000 SNe 
distributed uniformly in 0.1 < z < 1.7; this corresponds 
roughly (but not exactly) to what is expected from the 
SNAP space telescope [3a |. 

For a fixed dark energy model described by some values 
of the cosmological parameters wq and w a , let the prin- 
cipal components take values aj, with associated errors 
cr(aj) (note, we marginalize the results over the values of 
f?M) thus enlarging the errors in the PCs; in this way we 
can talk about models in the WQ-w a plane without further 
recourse to Om)- Further, from Eq. (fT0)l it follows that 
the principal component coefficients for the w{z) = — 1 
model are 
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where ei{z) is the shape of zth PC in redshift. We can 
now perform a simple x 2 test to determine whether w is 
constant: 
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where we have chosen to keep only the first M PCs, since 
the best-determined modes will contribute the most to 
the sum (note that the test is valid regardless of the value 
of M; in practice, our M is typically 3-5). Given that we 
have M degrees of freedom, each model (wq, w a ) will lead 
to x 2 which may or may be inconsistent with the ACDM 
model at some fixed confidence level. 

In order to determine the ability of future SN surveys 
to make this distinction, we focus here on the sensitiv- 
ity of the results to the measured magnitude uncertainty. 
We focus on this factor for two reason. First, we believe it 
will be the single most important determinant of the abil- 
ity of future surveys to possibly rule out a cosmological 
constant, and second, because our examination of past 
surveys suggests that one should consider the possibil- 
ity of redshfit-dependent magnitude measurement errors. 
This latter possibility is not unexpected. It is systemat- 
ically harder to determine the luminosity of ever fainter 
and more distant objects. 

Indeed, it is this possibility that originally motivated 
pfj| the current study. Measuring supernovae at ever 
higher redshift provide a useful lever-arm to distinguish 
between different equations of state, unless the magni- 
tude uncertainty increases more quickly with redshift 
than the redshift-distance relation diverges for differing 
values of the dark energy equation of state. For purposes 
of this example analysis, we used 217 SNla, from the 
initial large SN surveys [l|, Q , supplemented by several 
more recent measurements, fitting the quoted magnitude 
uncertainty as a function of redshift to a linear relation. 



Although the oft-claimed quoted magnitude uncertainty 
per low-redshift supernova is 0.15 we found that a free fit 
to the data tended to prefer a slightly larger value. Thus, 
for purposes of comparison, we considered both redshift 
independent uncertainties of 0.15 or 0.21 magnitudes, as 
well as redshift dependent uncertainties 0.15 + 0.079z or 
0.21 + 0.079z mag. 

Figure|3J shows the regions in the WQ—w a that are ruled 
out at 68% C.L. (region outside of the innermost closed 
curve) 95% C.L. (region outside of the middle curve) and 
99.7% C.L. (region outside of the outermost curve, coin- 
ciding with the blue shaded region). In other words, the 
blue region corresponds to dark energy scenarios where 
the ACDM model can be ruled out at > 3a confidence. 

In the four panels of Figure Q] we show cases when the 
error per SNa is 0.15 or 0.21 magnitudes, and when it is 
increasing with redshift as 0.15 + 0.079z or 0.21 + 0.079z 
mag. Clearly, increasing magnitude uncertainties with 
redshift can significantly affect the size of the indistin- 
guishable parts of model space. Note too that the shape 
of the region that cannot be distinguished from ACDM 
is of characteristic shape, is elongated roughly along the 
7(wo + 1) + w a rs direction — this is a direct conse- 
quence of the fact that the physically relevant quantity 

is w(z) — Wq + W a z/(1 + z). 

One other way to get a handle on whether w =/= —1 is 
to explore the most precise constraint on w that might 
be obtainable at any single redshift. For, if w =/= —1 at 
any redshift, this will establish unambiguously that the 
dark energy is not a cosmological constant. 

To explore the sensitivity of SNe for this purpose we 
parametrize dark energy equation of state as in Eq. Q , 
and study the accuracy in the equation of state w(z) at 
the best determined, or pivot, point. The pivot value is 
given by [Hil 



Cov(w , w a ) 

"pivot = W - — W a 

Cov(w a , W a ) 



(13) 



where Cov(cc, y) stands for elements of the covariance ma- 
trix element, while the pivot redshift at which w(z) is 
best determined is given by 



^pivot 



Cov(w , w a ) 



Cov(wo,w a ) + Cov(w a ,w a ) 



(14) 



Other parameter we are the matter density relative to 
critical, fij^ and the offset in the SNa la Hubble diagram, 
A4. Throughout we assume a flat universe. 

Using a Fisher matrix formalism, we estimate errors 
for a future SNAP-type survey [H] with 2800 SNe dis- 
tributed in redshift out to z — 1.7 as given by [36l |. 
and combined with 300 local supernovae uniformly dis- 
tributed in the z — 0.03 — 0.08 range. We study how the 
errors of the parameter of most interest, w p i vo t, change 
as the individual SNa errors vary. We again assume that 
the statistical error per SNa scales with redshift as 
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FIG. 4: Regions in the wo — w a that are ruled out at 68% CL. (region outside of the innermost closed curve) 95% CL. (region 
outside of the middle curve) and 99.7% CL. (region outside of the outermost curve, coinciding with the blue shaded region). In 
other words, the blue region corresponds to dark energy scenarios where the ACDM model can be ruled out at > 3<r confidence. 
This test is performed by "measuring" the principal components from each (wo,w a ) model and comparing them to those from 
the ACDM model. We show cases when the error per SNa is 0.15 or 0.21 magnitudes (top and bottom left panels), and when 
it is increasing with redshift as 0.15 + 0.079z or 0.21 + 0.079.Z mag (top and bottom right panels). 



a m {z) = a m (0) + Az, (15) 

where er m (0) = 0.10, 0.15 or 0.20, and we now let A vary 
in order to ascertain sensitivity to this parameter 

As can be seen in Fig. [5l the constraints on u> P i VO t are 
dependent on both a m (0) and A; we can fit the constraint 
via an approximate relation 

a(w pimt ) ^ 0.17 (2a m {0) + A). (16) 

This equation shows that the slope of the statistical 
error vs. redshift relation, A, contributes to cr(wpivot) one 
half as much as the intercept of the same relation, <r m (0). 



For example, a 20% increase in the SNa error between the 
redshift of and 1 (so that A = 0.2cr m (0)) will lead to a 
10% increase in error associated with the pivot value of 
the equation of state. 



IV. CONCLUSIONS 

SNe la are currently the strongest cosmological probes 
of dark energy, and are likely to remain the most solid 
source of information in the near future. This implies 
that control of the statistical errors is crucial. In par- 
ticular, if we are ever to address the central question of 
cosmology, 'Is the dark energy due to a cosmological con- 
stant?', we need to be able to unambiguously determine 
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FIG. 5: Constraints on the pivot value of the equation of 
state, uipivot, as a function of the intercept and slope of the 
magnitude error of SNe, <r m (0) and A. We have assumed a 
fiducial survey of 2800 SNe. 



that ro ^ -1 at some, or all redshifts. The extent to 
which SNe la measurements will allow such a determina- 
tion has been our prime concern in this analysis 

Figures 4 and 5 represent our primary results in this 
regard. Figure 4 makes clear that the model-independent 
constraints on w{z), for a two-parameter class of devia- 
tions around w(z) = — 1 (the two parameters are Wq and 
w a ] see Eq. ((3])), are quite sensitive to the measurement 
uncertainty in supernova magnitudes. Moreover, both 
Figure 4 and Fig. ^demonstrate that it is particularly im- 
portant to attempt to maintain control over measurement 
uncertainties for higher redshift supernovae. As Figure 4 



also demonstrates, allowing for possible variations in w 
with redshift implies to reduce the 95% confidence limit 
uncertainty in w(z) for a fit near wq = — 1 and w a = 
significantly below 10 % in wq will be challenging. 

This latter point is even more important when con- 
sidering whether one might utilize the measured value of 
w(z) at the pivot point to attempt to discern some time at 
which w =/= — 1. Redshift-dependent uncertainties in su- 
pernova magnitudes can easily almost double the inferred 
uncertainty in w at this point. Moreover, only if the 
planned large scale SNe surveys can maintain a uniform 
magnitude uncertainty per supernova, cr m , less than 0.10 
can we hope to derive a 95% confidence limit uncertainty 
in Wpivot of less than about 10% (see Fig. [5} which itself 
may not be sufficient to distinguish some non-standard 
dark energy models from a cosmological constant. 

It is not clear whether resolving the nature of dark 
energy is a 10-year or a 100-year problem. The answer 
partly depends on how much information measurements 
of dark energy can provide. Here we have addressed the 
former problem by studying the requirements of mag- 
nitude errors in a SNa la survey. Our results may be 
viewed as discouraging, or they may be viewed as in- 
spiring for those observers who enjoy a demanding chal- 
lenge. It is already becoming clear that, in order to have 
hope of detecting a deviation from the cosmological con- 
stant scenario, we need nature to be kind by producing 
a non-negligible deviation in w(z) at some redshift, and 
we need to control SNa magnitude errors, and systemat- 
ica in particular, to high precision. If either one of these 
requirements is not met, may have to rely on theorists 
to understand the nature of dark energy, which is option 
over which we may have much less control. 
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